Critical Phenomena in Head-on Collisions of Neutron Stars 
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We found type I critical collapses of compact objects modeled by a polytropic equation of state 
(EOS) with polytropic index F = 2 without the ultra-relativistic assumption. The object is formed 
in head-on collisions of neutron stars. Further we showed that the critical collapse can occur due to 
a change of the EOS, without fine tuning of initial data. This opens the possibility that a neutron 
star like compact object, not just those formed in a collision, may undergo a critical collapse in 
processes which slowly change the EOS, such as cooling. 
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Sec.l. Introduction and Motivation. In previ- 
ous studies [l|, Q we found that head-on colhsions of two 
neutron stars (NSs) could lead to prompt gravitational 
collapses, despite a conjecture to the contrary 0]. Upon 
collision, two 1.4 Mq neutron stars with a polytropic 
equation of state (EOS) could promptly collapse to form 
a black hole within a dynamical timescale, while NSs with 
smaller masses may merge to form a sing le NS ^. Fur- 
ther, in 0, we found that prompt collapses could occur 
even when the sum of the masses of the two colliding NSs 
is less than the maximum mass of a single stable equilib- 
rium star, in the case of NSs with the Lattimer-Swesty 
EOS. However, we were not able to determine the exact 
dividing line between the collapse and non-collapse cases, 
as a reasonably accurate determination would require nu- 
merical simulations with a resolution higher than what 
our 3D numerical code (GR- Astro could achieve on 
computers we have access to. 

In this work we investigate this dividing line. For the 
investigation, we have constructed an axisymmetric ver- 
sion of GR- Astro. GR-Astro-2D enables high resolution 
simulations of axisymmetric systems. Details of GR- 
Astro-2D and the numerical setup will be given in Sec. 2 
below. In Sec. 3, we report on the type I critical phenom- 
ena found at the dividing line for NSs with a polytropic 
EOS of F = 2. (Critical phenomena in gravitational col- 
lapse was first discovered by Choptuik [5|]; for review see 
e.g. [i,[3). In the super-critical regime, the merged ob- 
ject collapses promptly to form a black hole, even though 
its mass could be less than the maximum stable mass of 
one single NS in equilibrium with the same EOS. In the 
sub-critical regime, an oscillating NS is formed. Near 
the dividing line, the near critical solutions oscillate for 
a long time, before collapsing or becoming an oscillating 
NS. 

The exact critical solution presumably will oscillate 
forever, while neighboring configurations would eventu- 
ally deviate from the critical solution. In Sec. 4, we de- 
termine the critical index 7 as the time scale of growth of 
the unstable mode bringing a near critical solution away 



from the critical solution. We found 7 — lO.78(±O.6)Af0. 
This corresponds to a growth time of the unstable mode 
of about 0.054ms. In Sec. 5, we investigate the univer- 
sality of the phenomena with different critical parameter 
choices. Particularly interesting is the case in which the 
parameter is taken to be the polytropic index F, instead 
of a parameter associated with the initial configurations. 
The same critical index is found. In Sec. 6, we study 
gravitational waves emitted by the near critical solutions. 
We find that the gravitational radiations damp rapidly, 
suggesting that the unstable mode is not radiating. The 
determination of the possible existence of a non-spherical 
unstable mode/bifurcation point ([!]) requires a resolu- 
tion beyond the computational resources available to us. 

The existence of critical collapses in compact objects 
formed in collisions of NSs modeled by a commonly used 
EOS is by itself interesting. Further, the facts that 
(1) the growth timescale of the unstable mode is short 
(0.054ms), (2) it can appear in systems which deviate 
significantly from spherical, and (3) most importantly, a 
critical collapse could result through the softening of the 
EOS, suggest that critical collapses may actually occur 
in nature, e.g., a proto-NS initially supported by thermal 
pressure gradually cools down through neutrino emis- 
sions. When the EOS gradually changes in the timescale 
of seconds, the proto-NS evolves from the sub-critical 
regime toward the critical regime. The unstable mode, 
which grows in a timescale shorter than the evolution 
timescale of the EOS, may then be triggered, leading to 
a critical collapse. 

Sec. 2. GR-Astro-2D and the Numerical Setup. 

Simulations of NS collisions are carried out using the GR- 
Astro code, which solves the coupled system of the Ein- 
stein equations and the general relativistic hydrodynamic 
(GR-Hydro) equations as described in mm- The 
choice of numerical schemes as well as conventions used 
in this paper are given in [llj. A module GR-Astro-2D 
specially adapted to axisymmetric systems is constructed 
to provide the high resolution needed for this investiga- 
tion. GR-Astro-2D uses the same spacetime and GR- 
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FIG. 1; Lapse functions at center of 
collision vs time for NSs with slightly 
different masses. 



FIG. 2: Comparing the lapse functions 
and the 4-d scalar curvatures for the 
cases of lines 20 and 21 in Fig. [T] 



FIG. 
time. 



3: (a — Q:*)/a* v.s. coordinate 



Hydro evolution routines as GR- Astro but applies sym- 
metry conditions to restrict the simulation to a thin layer 
of spacetime (5 grid points across) containing the sym- 
metry axis . Initial data sets for head-on collisions of 
NSs are constructed by solving the Hamiltonian and mo- 
mentum constraints (HC and MCs) representing two NSs 
with equal mass separated by a distance along the z-axis 
and boosted towards one another at a prescribed speed. 
The numerical evolutions in this paper are carried out 
with the r freezing shift (or no shift) and the "1 -I- log" 
lapse (for details of the shift and lapse conditions, see 
[ll|). Extensive convergence tests have been performed 
including tests showing that the simulations are converg- 
ing throughout the evolution. In particular the critical 
index extracted (see below) has been shown to converge 
in a 1"* order manner, which is the designed rate of con- 
vergence due to our use of high resolution shock capturing 
treatment. 

NSs in this paper are described by a polytropic EOS: 
P = (F — l)pe with r = 2 (and cases close to 2). Here p is 
the proper rest mass density and e is the proper specific 
internal energy density. Notice that the "kinetic-energy- 
dominated" assumption has not been made, unlike ear- 
lier investigations of the critical collapses of perfect fluid 
systems (for review, see [1, Initial data sets are 

constructed with P = kp^ , where k = 0.0298c^/p„ (p„ 
is the nuclear density, approximately 2.3 x lO^"* g/crn^). 
For this EOS, the maximum stable NS configuration has 
an Arnowitt-Deser-Misner (ADM) mass of 1.46Af0 and 
a baryonic mass of I.6IM0. 

Sec. 3. Existence of the Critical Phenomena in 
Head-on Collisions of NSs. In the first set of sim- 
ulations, the two NSs are initially at a fixed distance 
(the maximum density points of the two NSs are sepa- 
rated by 3R, where R ~ 9.1Mq is the coordinate radius 
of the NSs). The initial velocities of the NSs are that 
of freely falling from infinity, determined by the Newto- 
nian formula plus the IPN correction. For the configu- 
rations investigated in Fig. [T] where the baryonic mass of 
each of the NS ranges from O.786M0 to O.793M0, the 
initial speed ranges from 0.15537 to 0.15584 (in units 
of c = 1). The computational grid has 323 x 5 x 323 



points, covering a computational domain of (vrr^ x height 
) = (tt X 38. 52 X 77.O)M0^ Each NS radius is resolved 
with 76 grid points, taking advantage of the octane- and 
axi-symmetry of the problem. 

Fig. H] shows the evolution of the lapse function a at 
the center of the collision as a function of the coordi- 
nate time, for systems with slightly different masses (all 
other parameters, including physical and numerical pa- 
rameters, are the same). The line labeled 1 in Fig. [T] 
(which dips to near t ~ 150Mq) represents the case 
of O.793M0. We see that after the collision, a promptly 
"collapses" to zero, signaling the formation of a black 
hole. Note that the total baryonic mass of the merged 
object I.59M0 is less than the maximum stable mass of 
a TOV solution of the same EOS in equilibrium. The 
prompt gravitational collapse of the merged object with 
such a mass indicates that it is in a state that is very 
different from being stationary [l|, Q . 

The fine labeled 41 in Fig.[T](which rises at t ~ I2OM0) 
represents the case where each of the NSs has the bary- 
onic mass O.786M0. The lapse at the collision center 
dips as the two stars merge, then rebounds. The merged 
object does not collapse to a black hole but instead form 
a stable NS in axisymmetric oscillations. The lapse at 
the center of the merged object oscillates around a value 
of 0.71, with a period of about I6OM0. 

For configurations with masses between the bottom 
line(l) and top line (41), the lapse a would rebound, 
dip etc., before eventually dipping to zero (a black 
hole is formed) or going back up (a NS is formed). 
The critical solution is found by fine tuning pc, the 
proper mass density as measured by an observer at rest 
with the fluid at the center of the star at the initial 
time. For the numerical setup used in the study, at 
around p^ = 6.128202618199 x 10""* (mass of each NS 
= O.79O7O949O26M0), a change of the pc by the 10*^ sig- 
nificant digit changes the dynamics from collapse to no 
collapse. In Fig. [1] we see that for these near critical con- 
figurations a oscillates at about 0.255 with a period of 
~ 40 A/0. As the lapse is given by the determinant of the 
3 metric, this represents an oscillation of the 3 geometry. 

For a more invariant measure, in Fig. [2] we plot as dot- 
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FIG. 4: Log{pc — Pc) vs. the departure 
time determined with Fig. [S] the slope 
gives the critical index. 
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FIG. 5: Lapses at the collision center 
for systems with the polytropic index 
r varying between 1.9997 and 2.0001. 



ted and long dashed lines the 4-D scalar curvature R at 
the collision center for two of the near critical solutions 
(lines 20 and 21 in Fig. [TJ they are the last ones to move 
away from the exact critical solution a,t t 3OOM0). We 
see that R oscillates with the same period as the determi- 
nant of the 3 metric (the lapse). As a collapses to zero, 
R blows up and in each such case we find an apparent 
horizon, indicating the formation of a black hole. Sim- 
ilar oscillatory behavior has been seen in other critical 
collapse studies d, 0, [1, 113] ■ We note that at late time R 
of the sub-critical case (line 21) tends to a small negative 
value as a static TOV star should. 

We note that while in Fig. [T]a change in the 10*'' sig- 
nificant digit of the total mass of the system can change 
the dynamics from that of sub-critical to supercritical, 
this does not imply that we have determined the criti- 
cal point to the 10*'' digit of accuracy. The exact value 
of the critical point is affected by the resolution of the 
numerical grid as well as the size of the computational 
domain. We have performed high resolution simulations 
with 76 grid points per R, (with computational domain 
covering 8.5i?), and large computational domain simula- 
tions covering 34i? (with resolution 38 grid point per R) . 
Convergence tests in both the directions of resolution and 
size of the computational domain suggest that the total 
mass of the critical solution in the headon collision case 
with the EOS given is at 1.58 ± 0.05Mq, with the error 
bound representing the truncation errors. 

Sec. 4. Determination of the Critical Index. 
The critical index 7 is determined through the relation 
T = 7log(p — p,), where T is the length of the coordinate 
time (which is asymptotically Minkowski) that a near 
critical solution with a parameter value p stays near the 
exact critical solution with p^. In Sec. 3 above, p is taken 
as the central density pc of the initial NSs. In Fig. [31 
we plot (a — a*) /a* at the center of collision against 
the coordinate time, where a* is the lapse of the critical 
solution to the best we can determine. Only the last part 
of the evolution is shown. We see explicitly the growth of 
the unstable mode driving the near critical solution away 
from the critical solution. We defined the "departure 
time" To. 05 as the coordinate time that a line in this 
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FIG. 6: The dominant 1—2, m=0 wave 
components. Line numbers here corre- 
spond to those in Fig. [T] 

figure reaches ±0.05 = ±5%. Likewise we define Tq.i, 
To. 15 and To.2- In Fig. 21 the departure times T0.05 and 
To. 2 are plotted against the log difference of p (taken to be 
Pc as in Fig. [1]) between the near critical and the critical 
solutions. With this, 70.05 defined as T'o.05/log(p — p*) is 
found to be 10.87, whereas 70.10, 7o.i5 and 70.2 are found 
to be 10.92, 10.93 and 10.92 respectively. We see that the 
value of the critical index does not depend sensitively on 
the definition of the departure point. 

Sec. 5. Universality. The above study uses the 
total mass/central density of the initial NSs as the critical 
parameter p. Next we fix the central density pc of the 
initial NSs at 6.12820305495 x lO^^'Afg-^ The initial 
coordinate separation between the center of the two NSs 
is fixed to be D = 27.5Mq. The initial velocity v is taken 
to be the parameter p. For each choice of v, the HC 
and MC equations are solved. Convergence with respect 
to spatial resolutions and outer boundary location has 
been verified. We find the same critical phenomena. The 
critical index is extracted in the same manner and found 
to be 10.78AfQ. . 

Other choices of parameter p have also been studied, 
including: (i) p ~ while fixing pc and v, and (ii) p = pc 
while fixing v and D. Note that the latter case is different 
from the case discussed in Sees. 3 and 4 above, where the 
initial velocity is determined by the free fall velocity up 
to the first PN correction. In all cases studied, we see 
the same critical phenomena with consistent values of 
the critical index 7. 

Next we ask: Is critical collapse possible only through 
fine tuning the initial data? If true, we would not ex- 
pect to see critical collapse phenomena in nature. We 
investigate the possibility of taking p = T, the adiabatic 
index, as slow changes of the EOS could occur in many 
astrophysical situations, e.g., accreting NSs and during 
cooling of proto-NSs generated in supernovae. We fix D, 
Pc V and vary F away from 2. The evolution of the lapse 
at the center of collision is shown in Fig. [5l We see be- 
havior similar to that of Fig. [H The critical index 7 is 
found to be again 1O.78M0, consistent with the values 
found by fine tuning the initial configurations. 
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Sec. 6. Gravitational Wave Signals Prom Near 
Critical Collapses. Can near critical collapses have a 
signature in their gravitational radiations that one can 
identify in observation, in view of the possibility that 
near critical collapse of a neutron star like object may 
occur in nature? 

In Fig. [6l we show the dominant gravitational wave 
component (even parity L = 2,m = component of the 
Moncrief gauge invariant Q [l5| , z-axis is the symmetry 
axis) for four cases shown in Fig.[T] For the line 41, the 
merger leads to a NS oscillating in a non-spherical fashion 
and emits gravitational waves. For the line 1 the merged 
object moves away from the critical solution at around 
t ~ 130 and collapses to a black hole. The asymmet- 
ric collapse generates gravitational wave. The simulation 
was terminated at i ~ 330 as the metric functions can no 
longer be adequately resolved with the formation of the 
black hole. In contrast, the radiations from the near crit- 
ical solutions represented by lines 20 and 21 (the dotted- 
dashed line and the solid line respectively, nearly coincid- 
ing with one another all the way) decrease in amplitude. 
At around t ^ 270, the unstable mode sets in and the 
two near critical solutions 20 and 21 move away from the 
critical solution (as shown in Fig. [1]). However, there is 
only a negligible amount of wave emitted, independent 
of whether the unstable mode leads to a black hole (line 
20) or an oscillating NS (line 21). This suggests that the 
unstable mode may be a spherical mode that does not ra- 
diate. The waveforms reported here suffer from that the 
computational domain (covering only ~ 1 wavelength) 
cannot be extended further out, limited by the compu- 
tational resources available to us. Further investigations 
of the waveforms will be reported elsewhere. 

Sec. 7. Conclusion. We found that the dividing 
line between prompt and delayed collapses of a compact 
object formed in a head-on collision of NSs with a poly- 
tropic EOS occurs at a mass less than the maximum sta- 
ble mass of a single equilibrium star. The EOS is one 
that is often used in modeling NSs. There exists a type I 
critical collapse phenomena at the dividing line. Univer- 
sality of the phenomena with respect to different choices 
of the critical parameter is confirmed and the critical in- 
dex extracted. The growth time of the unstable mode 
which brings a near critical solution away from the crit- 
ical one is found to be 0.054tos, for the EOS studied. 



The study suggests that, despite the highly asymmetric 
initial data, the final critical collapse may not generate 
gravitational radiation. 

We found that critical collapses could happen without 
fine tuning of initial data, but instead, through a gradual 
change of the EOS. This opens the intriguing possibility 
that, e.g., when a proto NS cools and loses thermal sup- 
port on a timescale longer than 0.054ms, the collapse 
could exhibit critical behavior. This and related ques- 
tions will be investigated elsewhere. 
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